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Abstract 

The rate for the production of a pair of massive fermions in e + e~ annihilation plus real or virtual 
■ radiation of a pair of massless fermions is calculated analytically. The contributions for real and 

virtual radiation are displayed separately. The asymptotic behaviour close to threshold and for high 
energies is given in a compact form. These approximations provide arguments for the appropriate 
choice of the scale in the 0(a) result, such that no large logarithms remain in the final answer. 

o\ : 

►> ', 1. Introduction 

\ QED reactions leading to four fermion final states /1/1/2/2 with fermion masses nil and 777-2 have 
been considered in the literature a long time ago |T| and the final result was at that time expressed 

in 1 
o 

context of QED. However, it is fairly evident that a large part of the results can easily be transferred 
to the corresponding QCD reactions. The analogy becomes even closer when considering the mass 



in the form of a two dimensional integral. The cross sections were calculated for leptons in the 



assignments of interest for actual physical reactions: either equal masses mi = 777,2 or alternatively 
mi 3* 777,2- This mass hierarchy applies equally well to leptons and to quarks and hence to QCD 



Oh 
1 

qj , calculations and will be exploited in this work. 

Apart of the "exclusive" channel, where all four fermions are detected, also the inclusive rate 
for fxfi + anything is of practical interest. To 0{a) this calculation has been presented in the classic 
^ ■ book by Schwinger [Q]. Higher order results are available in the context of QCD in the limit of 
vanishing fermion mass [||, |]] or including mass corrections through an expansion in m 2 /s. Terms of 
order m 2 /s where calculated up to ||, the m 4 /s 2 terms up to a 2 , ||. This expansion is, however, 
inadequate close to threshold. In this kinematical region the full calculation up to order a 2 would be 
required. This is particularly desirable in view of the fact that the ambiguity in the scale u 2 in a s 
leads to a large uncertainty in the leading order correction. Close to threshold both the mass and the 
three- momentum of the produced fermions seem to be reasonable choices for u 2 , giving rise, however, 
to drastically different predictions. 

The analytical calculation of the cross section, including the full mass and energy dependence 
and all (real and virtual) gauge boson contributions seems like a difficult task. However, the subclass 
of diagrams involving the real and virtual radiation of light fermions is more accessible. In this 
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paper the result for both real and virtual radiation will be presented for arbitrary mf and s, in the 
limit m\ S> m?,. Combining the two contributions, one arrives at a result which still exhibits mass 
singularities of the form h^ml/s). They can be removed by adopting the MS scheme. Normalizing 
the coupling constant at scale fj? = s eliminates all large logarithms, at least away from the threshold 
region. This provides the first step in the calculation of corrections to order a 2 . The relatively compact 
analytical result can then be studied in the limit close to threshold as well as in the high energy region. 



2. Real radiation 

For defmiteness the cross section for the production of the /1/1/2/2 final state in e + e~ annihilation 
through a virtual photon will be considered, normalized relatively to the point cross section. The 
result is evidently also applicable to Z decays through the vector current. The relevant Feynman 




Figure 1: Characteristic Feynman diagrams describing the production of a pair of massive and a (real 
or virtual) pair of massless fermions. 



amplitudes can be derived from the four fermion cuts of the diagrams with two closed fermion loops, 
with representative examples depicted in Fig. |l[ The rate can be expressed by a two dimensional 
integral M: 
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where 



A(l,y,z) = l + y 2 + z 2 -2{y + z + yz). 



(2) 



(3) 



For arbitrary mi, rri2 and s already the first integration leads to elliptic functions and fairly lengthy 
expressions. However, in the limit m\ 3> m| the radiation can be split into two parts: "soft" radiation 
with energy of the ^2/2 system smaller than a cutoff A, with m2 <C A <C mi, and the remainder, 
denoted "hard" radiation. The two parts can be integrated separately, and their sum is given by 



^ = /^ln 2 ^ + /«lu^ + /W 



(4) 
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As expected, the cutoff A cancels in the sum. The functions multiplying the second and first power 
of the logarithm are closely related to the well-known Schwinger result [Q] for real radiation of a light 
vector particle with mass A <C mj: 
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The evaluation of the function without logarithmic enhancement constitutes the main effort of this 
section. The three functions f^' 1 ' 2 ^ are given by 
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Li 2 , Li 3 denote the di- and trilogarithms, £(2) and ("(3) the Zeta function of the respective argu- 
ments |0]. The details of this calculation will be given elsewhere g. This result is directly applicable 
for example to the annihilation of e + e~ into t + t~ with the additional "final state" radiation of e + e~ 
(or /i + p~) which adds incoherently to the process with r radiated from the light fermion final state ||. 
The behaviour close to threshold (w — ► 0) and in the high energy region (m 2 /s — ► 0) is of relevance 
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for direct applications and for the discussion of our results presented below. One finds 
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Close to threshold radiation of light fermions is strongly suppressed, similar to the radiation of photons 
as calculated in O{o). The leading logarithms in the high energy limit coincide with those given in ]!]]. 
The prediction for Rf 1 f 1 f 2 f 2 (with toi/to-2 = m T /m e = 3477.5 chosen for illustrative purpose) based 
on eq. (|l]), is shown in Fig. ^ (solid line). Also shown are the high energy approximation, eq. (11), 
(dashed line) including the linear term in x and the threshold approximation, eq. (|To|) , (dashed-dotted 
line). It is evident that the high energy approximation can be used for values of x = mf/s between 
and about 0.125 corresponding to values of w from 1 down to 0.7 and hence surprisingly close to the 
threshold. For w below this value the threshold approximation provides an adequate description. 
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3. Virtual corrections 

Virtual corrections in the present context arise from the two particle cut of the "double bubble" 
diagrams (Fig. [J). The 0(a 2 ) corrections SA^ to the lowest order vertex can be classified into contri- 
butions to the Dirac (Fi) and the Pauli form factor (F2) 
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where a^ v = 5(7^,71/] and g denotes the photon momentum flowing into the vertex. Using the 
dispersive methods applied already in [|l0|, [yj the calculation of can be easily reduced to a one 
dimensional integration. It results from the convolution of the massive vector boson exchange vertex 
correction to f\f\ production with the absorptive part of the vacuum polarization of fermions with 
mass m2. Denoting the vector boson mass by A the convolution reads 
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The QED normalization -Fi(O) = is evidently adopted. The evaluation of eq. (Op is tedious for 
arbitrary m\ and m 2 and will be presented in detail in ||]. In the limiting case m\ 3> mi, however, 
the result is drastically simplified: 
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Again the functions /{ and are closely related to the very well-known logarithmically divergent 
and constant pieces of the corresponding one loop corrections for small A 0]: 
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The similarity of these relations with eq. (|5|) is evident. Interesting special cases are again the behaviour 
close to threshold and for high energies. Let us discuss the former: 
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The Coulombic behaviour ~ 1/w is evident from this result. For F\ the Coulomb singularity is 
modified by the logarithmic factor ln(m2/s) which is responsible for the "running" of the coupling 
constant in the 0(a) result. It is instructive to combine the vertex correction of O{o) and 0(a 2 ) in 
the region close to threshold. As an illustrative example we will examine the Dirac formfactor for this 
case. The infrared divergent part of F\ vanishes for w — > and therefore 
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The fine structure constant a, defined at vanishing momentum transfer, is related to the MS coupling 
constant at subtraction point /i 2 by 
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At this point it becomes obvious that the natural scale for a^g in the threshold region is given by 
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as far as the 1/w terms are concerned. [This holds true as 



long as w ~ a. Below this value the approximations used in this work are no longer applicable.] For 
the correction resulting from transverse photon exchange, which are not enhanced by 1/w, the scale 
fj? = s/4 is appropriate. This suggests the following form of the Dirac formfactor in the threshold 
region 
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where the scale in the 0(a 2 ) term is not yet determined. It is clear that the virtual corrections will 
dominate the rate close to the threshold. 
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For high energies, on the other hand, one finds 
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For the special case mi = m 2 = m the integrals dl4] ) lead to particularly simple results 
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The logarithmically enhanced and the constant parts are in agreement with 

Finally, the 0(a 2 ) contribution of the virtual corrections to the rate for the case mi 3> m 2 is 
given by 
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4- The total rate 

Combining real and virtual radiation one thus arrives at 

9 R + 9 V = -\w\a^ + ^ +W (3- W 2 ){f?Hf^) +W "ff ) . (40) 

The quadratic logarithm in m 2 /s from the real and the virtual radiation cancel. A linear logarithm, 
however, remains. Its origin can be easily understood through the running of the coupling constant a. 
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The prefactor W is identical to the correction function of O(a) derived by Schwinger [Q]. Therefore 
the expression for the rate for the inclusive /i/i final state including O(a) photonic corrections plus 
photonic 0{ot 2 ) corrections due to one light fermion with mass m,2 reads 
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[Note that the massive quark is not accounted for, consistent with the fact that virtual heavy fermion 
loops are not considered in eq. (Q). Adding virtual corrections eq. ( pioj , |3G| ) and the real radiation, 
e.g. based on a numerical evaluation of eq. (|l|) one would thus include "double bubble" diagrams with 
two massive fermions. This will be done in Q.] 

Relating again the fine structure constant a to the MS coupling at the scale p 2 , the mass 
singularities disappear and one finds 
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The behaviour close to threshold for the choice p 2 = s is easily read off from eq. (42): 
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The discussion following eq. (|30| ) applies equally well to this formula, since real radiation vanishes 
close to threshold. 

In the high energy region one finds 
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In Fig. ^ the comparison between the exact 0(a 2 ) correction for u 2 = s (solid line), threshold ap- 
proximations (dashed dotted lines) and high energy expansions (dashed lines) is performed. Eq. ( |44"D 
provides an important consistency check on our result. It is straightforward to relate pole and MS 
definition for the remaining fermion mass mi taking again into account in 0(a 2 ) only the contribution 
from one light virtual fermion: 
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Replacing the pole mass by the running mass at the scale fj, 2 = s the logarithmic factor of the m 2 
term disappears as expected from general considerations. The structure of the logarithms of the m 2 



and m 4 terms coincides with the expectations from 
by the proper SU(3)-coefficients one obtains 



In fact, after replacing the abelian factors 
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where now the number of light fermions, nj, is displayed explicitly. The relation to the rif dependent 
terms of eq. (27) in || is evident. 
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Figure 3: 0{a 2 (s)) correction to the inclusive production rate Rf 1 f 1 based on the exact result (solid 
line) and approximations described in the text. 
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5. Summary 

The rate for the production of a pair of massive fermions in e + e~ annihilation plus real and virtual 
radiation of a pair of light fermions has been calculated analytically. This result, together with || 
can be considered as a first step towards the evaluation of the production cross section for heavy 
fermions in 0(a 2 ). The expansion of the result for energies close to threshold and for high energies 
and subsequent comparisons with earlier asymptotic formulas provide important cross checks. The 
transition to the MS scheme leads to interesting insights into the proper scale of the coupling constant. 

Acknowledgement: We would like to thank K. Chetyrkin and M. Jezabek for helpful discussions. 
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